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Abstract: We study the effect of radiative corrections on the structure of neutrino mass 
matrix. We analyze the renormahzation of the matrix from the electroweak scale niz to 
the scale mo at which the effective operator that gives masses to neutrinos is generated. 
Apart from Standard Model and MSSM, non-standard extensions of SM are considered 
at a scale mx intermediate between mz and tjiq. We find that the dominant structure 
of the neutrino mass matrix does not change. SM and MSSM corrections produce small 
(few percents) independent renormahzation of each matrix element. Non-standard (flavor 
changing) corrections can modify strongly small (sub-dominant) matrix elements, which are 
important for the low energy phenomenology. In particular, we show that all sub-dominant 
elements can have purely radiative origin, being zero at mo- The set of non-zero elements 
at mo can be formed by (i) diagonal elements (unit matrix); (ii) Mee and M^r] (iii)Mee and 
//r-block elements; (iv) //r-block elements. In the case of unit matrix, both atmospheric 
and solar mixing angles and mass squared differences are generated radiatively. 
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1. Introduction 

Experiments on neutrino oscillations Q, neutrinoless 2(3 decay and tritium /5 decay Q 
give information on neutrino masses and mixing angles. In principle, also CP violating 
phases can be measured. Using this information, it is possible to reconstruct (at least 
partially) the Majorana mass matrix of neutrinos at the electroweak scale mz- 

The origin of the neutrino mass matrix is, most probably, in new physics above a very 
high scale mg, where the lepton number is violated. To find the structure of the matrix at 
mo, it is necessary to take into account the renormalization effects between mz and mo. 

The Renormalization Group Equation (RGE) for the neutrino mass matrix has been 
extensively studied both in the Standard Model (SM) and Minimal Supersymmetric Stan- 
dard Model (MSSM) The RGE has been also considered, in the context of the 
see-saw mechanism, for non-degenerate heavy right-handed neutrinos |l^. Low energy 
threshold corrections have been computed; their effect can be important in the MSSM 
case 1 11]. A non-standard source of neutrino masses (operators in the Kahler potential) 



has been considered in [12| and the corresponding radiative corrections have been studied. 
The goal of these analyses was to understand how radiative corrections could affect the 
mass squared differences Am?^ and the mixing angles 9ij observed in neutrino oscillation 
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experiments. A general conclusion is that the effect of running is small for hierarchical 
mass spectra, while, in the case of quasi-degenerate mass spectrum, the observables can be 
strongly modified by radiative corrections. In particular, in the degenerate case Am^^j can 
have purely radiative origin . 

Mass squared differences and mixing angles are the outcome of the diagonalization of 
the mass matrix and it is this matrix the object more closely related to the underlying 
theory. In contrast with previous studies, we will analyze in detail renormalization effects 
on the mass matrix structure. 

In section § we discuss the generation of neutrino mass matrix at the scale mo. We 
study first the radiative corrections in SM and MSSM (section ^). Then, we consider other 
extensions of the SM at some intermediate scale mx, with mz < mx < niQ (section We 
identify the features of the mass matrix which can be explained by radiative corrections. 
Finally (section |5|) , we consider four specific matrix structures at ttiq: (i) the matrix propor- 
tional to the unit; (ii) the matrix with only ee and /it elements different from zero; (iii) the 
matrix with non-zero M^e and /ur-block elements; (iv) the matrix with non-zero //r-block. 
We will study the predictions for the low energy parameters in the case of standard and 
non-standard radiative corrections. 



2. Generation of the neutrino mass matrix and scale mo 

With the SM fields, one can construct a unique Lorentz and gauge invariant effective five 
dimensional operator that gives Majorana masses to neutrinos [14|: 

^(L^a ia2 0)(0icT2L^) + h.c. (2.1) 
mo 

Here Ca/s are dimensionless couplings, mo is the mass scale at which this effective interac- 
tion is generated, a, (3 are flavor indexes, (j) and L are the SM Higgs and lepton doublets. 
An analogue operator appears in the case of MSSM (without R-parity violation). After the 
electroweak symmetry breaking, the operator (|2.l| ) generates the neutrino mass matrix: 



M, 



al3 



mo 



Different possible mechanisms can lead to the effective operator ( |2.1[ ). One possibility 
is the exchange of heavy right-handed neutrinos (type-I see-saw [^). In this case mo 
should be identified with the mass of the lightest right-handed neutrino m/j and one gets 

C^f, = -^{Y,M^'YJU, (2.2) 

where is the matrix of neutrino Dirac Yukawa couplings and Mpi is the Majorana mass 
matrix of right-handed neutrinos, both evaluated at the scale mo = niR. 

Another possibility is to introduce a scalar isotriplet A with hypercharge 1 and renor- 
malizable coupling to the SM lepton doublets: 

Y^f^l7c,a'LpA, + h.c. . (2.3) 
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A heavy triplet {uia 3> mz) can get a small induced VEV due to the coupling M/^fpA*(p(p+ 
h.c. (type-II see-saw [16|): 



If the exchange of the triplet is the dominant contribution to the neutrino mass matrix, 
tuq should be identified with rriA and 

Cc.0 = Y^^^. (2.4) 
mo 

The following remarks are in order: 

1. The neutrino mass matrix can receive both contributions of the form (|2.2| ) and (p.4|). 
In this case mo = min{m/{, rriA}- 

2. The running between mo = min{m/j, m^} and mo = max{m/j,mA} can be impor- 
tant for the neutrino mass matrix structure, if Nji and A contributions are compa- 
rable. 

3. There can be a hierarchy among the masses of right-handed neutrinos, so that the 
operator ( |2.1| ) is formed in a large energy interval. Corrections in this interval can 
be important |jlO|. 



The relative size and the structure of the various {Nji, A, maybe some other) contributions 
to the neutrino mass matrix are model-dependent and deserve a separate analysis. In this 
paper we will study the running below mo- 

3. Renormalization of the mass matrix in the SM (MSSM) 

Let us consider, first, the renormalization of the mass matrix when the only particles 
with mass below mo are the SM (MSSM) particles. The /?-function of the operator (p.lj), 
Pm = /^ajl-^! be written as ^ |5|: 



+ KsmM , 



167r^/3; 



2 oMSSM 



M 



M{Y^Yi) + [Y^Yifu] + KmssmM , (3.1) 



where Yi is the 3x3 matrix of charged lepton Yukawa couplings and ^i^5M{M5S'M) is a real 
parameter describing flavor universal radiative corrections. 

The flavor non-universal corrections (terms in square brackets in eq. ( |3.lD ) come from 
two types of diagrams only (figure |^), generated by charged lepton Yukawa couplings: 
renormalization of the wavefunction of lepton doublets (a) and vertex correction (6). In 



the SM, the coefficient —3/2 in eq. (3^) is the sum of a contribution 1/2 from (a) and —2 
from (b). In the MSSM, only (a) contributes, due to SUSY non-renormalization, and the 
coefficient is 1/2 x 2 = 1, where the factor 2 corresponds to the double number of particles 
with respect to SM. 
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(a) 



(b) 





L * L 

Figure 1: The flavor non-universal diagrams in the case of SM (MSSM). 



In flavor basis, the matrix Yi is diagonal and it can be made real by reabsorbing the 
phases in the fields: 

Y.^Yi = dmg{yl yl,y',). 

Furthermore, ya, taken real at a certain scale, will remain real during the running [^. In 
this basis, the RGB equation for M^^ can be easily integrated, giving Q: 



where 



tz 



^ log(mo/m^) 
167r2 



-k f\yl{t) + yl{t))dt 
Jo 



Mai3{mo) , 



(3.2) 



l (SM) 



k = 1 (MSSM) . 



The flavor universal corrections are contained in the prefactor 

rtz 



Ik = exp 



K{t)dt 







Ik does not change the structure of the mass matrix and, moreover, the overall renor- 
malization effect is small: ~ 1. We will consider only 1-loop radiative corrections and 
therefore neglect the evolution of ya in eq. ( |3.2D , taking the values of charged lepton Yukawa 
couplings at the electroweak scale: ya{t) ~ ya(0). 

Several important conclusions follow immediately from eq. ( |3.2D : 



• In flavor basis, each element of the mass matrix evolves independently from the values 
of other elements. The value of the element at mz is proportional to its value at mg: 

Map{mz) oc Mapimo) . (3.3) 

In particular, elements which are zero at some scale, remain zero at any scale. 



The phases of M^p do not evolve, because the r.h.s. in eq. ( |3.2D is real. In contrast, 
the three physical CP violating phases depend on the moduli of M^p and therefore 
can be strongly renormalized. However, if there is no CP violation at some scale (M 
is real), no CP violation can be induced by radiative corrections at any scale. 
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• The corrections to matrix elements have opposite sign in SM and MSSM, due to the 
different sign of k. The same matrix at mz wiU develop different features at the scale 
uiQ depending on whether low energy supersymmetry is present or not. 

Notice that, in principle, SM and MSSM could be discriminated by the ordering of mass 
eigenvalues: exactly degenerate neutrinos at the scale mo could be split into a normal 
mass spectrum in SM and into an inverted spectrum in MSSM or vice versa, depending on 
mixings and phases. However, using only SM (MSSM) radiative corrections, it is hard to 
reproduce data starting with exactly degenerate neutrinos [0] . Also the scenario with only 
two degenerate neutrinos at high energy has been recently studied [^] and the difference 
of predictions between SM and MSSM has been analyzed. 

The largest flavor dependent contribution to the running in eq. (|3.2[) is due to the r 
Yukawa coupling: 

^ « 10-2 (SM) 

, (3.4) 

" tan/3- 10-2 (MSSM) 



yAmz) 



V cos (3 

where v ~ 246 GeV is the VEV of the SM Higgs and tan /? is the ratio of VEV's of the two 
MSSM Higgs doublets. The other Yukawa couplings are much smaller, therefore the largest 
corrections are for elements which have the r-flavor. Neglecting ye and corrections, 
eq. (^.21) gives: 

Mrr{mz) ~ lKMrr{mo){l - 2/ce^) , 

Merimz) ~ /E-Mer(mo)(l - kCr) , 

M^r{mz) ~ /i^M^^(mo)(l - /ce^) , 
where, taking mo = 10" GeV, 

_ vlimz) , rriQ . 2 /o r:\ 

log < Aan/3\ /^.ooTv.N • (3-5) 




167r2 ^ mz « 3.5(n - 2)10-^ (MSSM) 



The effect of running is apparently very small even for M-j-r- 

We conclude that radiative corrections (in SM and MSSM) practically do not change 
the structure of neutrino mass matrix up to the scale niQ. Zero elements remain zero and 
non-zero elements acquire very small relative corrections, which are at most few percents. 
Symmetry properties of the matrix at mo are almost unchanged by running to low energy, 
where the matrix structure can be reconstructed using experimental input for mixing angles 
and mass eigenvalues, as well as CP violating phases |]l^, In other words, the struc- 
ture of the mass matrix is stable with respect to radiative corrections in SM and MSSM, 
independently of type of mass ordering (normal or inverted), level of degeneracy, values of 
mixing angles and CP violating phases. 

3.1 Radiative corrections and observables 

Even though the matrix structure is stable, observables, i.e. the values of the mass dif- 
ferences and the form of the mixing matrix, can be strongly affected by the radiative 
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corrections. There is a number of studies in which the RGE's for mixing angles 

and Am?j have been analyzed and conditions for strong renormalization effects have been 
identified. 

In contrast with previous studies, here we discuss the effect of renormalization on 
observables in terms of mass matrix. The matrix of corrections can be written as: 

AM^-eKM^'-k ... 2e^M0^ (e^ + e^)^^^ , (3.6) 

where M° is the matrix at the scale mo and ^K,e,^i can be obtained from the expression (^.5|), 
by substituting with K,y'^,y'^ respectively. 

The effect on a given observable depends on how strong is the influence ( "imprint" ) of 
this observable on the structure of the matrix. When two eigenstates are almost degenerate 
in mass, the matrix structure depends very weakly on their mass squared difference Amjj 
and on the mixing angle 6ij []T9| , [20[| . In other words, Amfj and Oij are not "imprinted" in 
the matrix structure. In this case the small radiative corrections ( |3.6D can strongly modify 
these observables. Let us give a rough estimation: 9ij and Amf^ can receive large radiative 
corrections if the absolute neutrino mass scale m is large enough to satisfy Am?^ ~ e,- m?. 
Because of the smallness of Cr, this condition requires quasi-degenerate neutrino masses. 

In section |5|, we will use eq. ( |3.6D to discuss the effect of radiative corrections for some 
particular structures of M^. 

4. Renormalization of the mass matrix due to new particles 

It may happen that new physics exists at some intermediate scale mx in the range mz—rriQ, 
which does not contribute to neutrino masses at tree-level but leads to renormalization 
effects. 

4.1 Non-standard Yukawa interactions 

Let us consider new fermions and scalar bosons, with mass ~ rnx- We assume that 
new scalars have zero VEV's, otherwise they would generate neutrino masses at the scale 



mx < TnQ. We analyze radiative corrections to the operator (2.1) induced by the couplings 
of these scalars and fermions to the lepton doublets La- 

• An extra scalar doublet, cj)', with the coupling 

+ h.c. , (4.1) 

contributes to the wavefunction renormalization (diagram analogue to figure [l| 



There is no vertex diagram with (p' in the loop (analogue to figure |l^), because only 
(j) enters the operator (p.l|). 



An extra scalar singlet p or triplet x can couple with two lepton doublets: 

Y^^Laia2L''pP + h.c., (4.2) 
Y^^Laia2a'L'pXi + h-c . (4.3) 
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(a) 



(b) 





Figure 2: The diagrams generated by a new scalar singlet p or triplet x (eqs. (4.2), (4.3)). 

P' X 




Figure 3: The diagram generated by a new scalar singlet p (triplet x) and a new fermion doublet 



Dr (eqs. (4.4), (4.5)) 



The renormalization of the operator ( |2.1D is induced by the two diagrams in figure |2|. 
In the singlet case, the matrix Y is antisymmetric, because Laicr2L'^ = —Lpia2L'^ 
(in the triplet case, Y is symmetric). 

If a new right-handed doublet fermion D^i exists at the scale mx, p and x can have 
the following interactions: 



Y^PLo^DnP + h.c, 
Y^^L^aWnXi + h.c. 



(4.4) 
(4.5) 



In this case there are no vertex corrections, because D_r and p (x) do not enter in 
the operator (p.l|). Only the diagram shown in figure contributes. 



A right-handed fermion singlet Sr, with the Yukawa interaction 



Y^LaCpSn + h.c. . 



(4.6) 



contributes both to wavefunction and vertex renormalization, as shown in figure |^. 
The scalar doublet in the diagram ^(a) can be the SM one, (p, or some new doublet (/>'. 

The new chiral fermions {Dji, Sn) can lead to anomalies, however one can consider them 
as part of new vector- like fermions. The vertex diagrams in figures ^ and ^ give no 
contributions in the supersymmetric case, because of non-renormalization theorem. 
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For any of the interactions in eqs. ( [4.1[ )"( |4rq ), the contribution to the /3-function of 
neutrino mass matrix elements can be written as 



I67:\pl^)af3 = k^Y^ M{Y^Y) + {Y^YfM +4'^ YM^Y^ . (4.7) 



.(2) 



Q/3 



(i) 

where Y are the new particle Yukawa couplings and the prefactors ky depend on the 



type of particles and interactions considered. The first term in square brackets of eq. 
(analogue to those in eq. (^)) corresponds to all diagrams discussed above (analogue to 
those in figure |T|), except the diagram in figure |2|6. This diagram generates the second 
term in square brackets of eq. ( |4.7| ). Therefore, ky^ ^ only for the interactions in 
eqs. O, (13). 

Let us emphasize the main differences in the running with respect to SM (MSSM): 



If the matrix Y^Y (or, in the case of eqs. ([4. 2D, (4.3), the matrix Y) is not diagonal 



in flavor basis, the RGB's for different matrix elements are coupled. A given matrix 
element receives corrections proportional to other matrix elements; in this way small 
elements can be modified significantly. 

The size of corrections depends on the size of Yukawa couplings Y. It can be much 
larger than in SM. In the perturbative regime, y < 1, the effect can be of the order 
of few percents, as in MSSM. 

The corrections can be further enhanced if several non-standard multiplets are 
present. E.g., one can introduce three generations of new particles. In this case 
corrections can be as large as ~ 10% of the large matrix elements. 

The size of corrections depends on mx- Corrections are suppressed for mx ~ mo- 

In all diagrams but figure ^6, only one external lepton leg is involved in flavor changing 
interactions. As a consequence, a given element Ma/s receives contributions propor- 
tional to matrix elements in the rows (columns) a and /3 only. In contrast, if the 
diagram in figure ^ is present, all matrix elements can contribute to the renormal- 
ization of a given element M^^. 
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Since radiative effects are small, with a good approximation we can consider only lowest 
order corrections to matrix elements. Using eq. ( [4.7D , we get: 



M{Y^Y) + {Y^Y^m] + k^Y^ \yM^Y^] ) . (4.8) 

J a/3 L J a/3 J 



4.2 Non-universal U(l) gauge interaction 

Let us consider the effect of extra heavy vector bosons X^, which have flavor non-universal 
interactions. These bosons can be related to the existence of horizontal (flavor) gauge 
symmetries [pT| . 

We restrict ourselves to the case of an extra U(l)x gauge group. In the flavor basis, 
the interaction of lepton doublets with the new gauge bosons has the form 

gxQa/3Lal^X^Lf3 , (4.9) 

where gx is the gauge coupling and Q is the hermitian matrix of "charges" , which can be 
non-diagonal in flavor basis. The possible anomalies of the extra U(l)x gauge group can 
be canceled using the Green-Schwarz mechanism [^]. 

In figure ^ we show the two 1-loop gauge diagrams that give non-universal radiative 
corrections to the operator ( |2.1| ). Their contribution to the /3-function of matrix elements 
takes the following form: 

16vr2(/3f,)„^ = g^xkPlM + {Q^)'' + g\k%\Q^ M QU . (4.10) 

The first and the second square brackets of eq. ( |4.10D are the contribution of the diagrams 
in figure and |5|J), respectively. The /3-function ( 4.10| ) has analogous features to the one 



in eq. ( [4.7[ ): the RGE's of different elements are coupled because of the off-diagonal entries 
in the matrix Q; the corrections are proportional to the gauge coupling g\ and cannot 
be larger than few percents of the largest element in M. Like in figure ^6, in figure ^6 
both lepton external legs enter in the loop. As a consequence, a given element Ma/3 can 
receive contributions proportional to all matrix elements. In lowest order, corrections can 
be written as: 

^ _ gilog(moMx) ^^ia)^^ ^ {QYMU + kflQ^M QU) . (4.11) 
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5. Radiative generation of the sub-dominant matrix elements 

The values of mass squared differences and mixing angles, measured in neutrino oscillation 
experiments are (90% C.L.): 

Am^^i = Amf^ = (7 tf) ■ IQ-^eV^ , tan^ Ou = 0.42+°;^ (LMA MSW) ; 

Aml^ = Amig = (2.5 t^) ' W~^eY^ , tan^as = 1 t^oH ! sin^ig < 0.2. (5.1) 

Using these data, one finds that the neutrino mass matrix in flavor basis can have a 
hierarchical structure, with some elements much smaller than the others. Small elements 
are suppressed by factors S13, Am^^J Am'^^^ or, in the case of degenerate spectrum (mi « 
1712 ^ "^3)) by Am^tjjj/m^. All possible hierarchical structures allowed by the data have 
been identified in p0[| . 

Non-standard radiative effects can generate non-zero matrix elements even if they are 
zero at the scale niQ. These elements can receive a radiative contribution up to ~ 10% of 
the largest matrix element. In what follows we will consider mass matrices with various 
hierarchical structures and study the possibility to generate all small elements of these 
matrices radiatively. It has been found that matrices with three or more exactly zero 
elements at the scale mz cannot reproduce the experimental data [p^ j. We will show 
that these matrices, realized at mo, agree with phenomenology if non-standard radiative 
corrections are taken into account. 

Let us write the neutrino mass matrix with hierarchical structure as 

M = Md + Ms , 

where Md is the matrix of dominant elements and M5 is the matrix of sub-dominant 
elements. In general, large matrix elements are the sum of a dominant contribution from 
Mjy and a small correction from M5. Small matrix elements are contained in M5 only. We 
assume that, at the high scale mo, 

M(mo) = MD(mo) , Ms(mo) = (5.2) 

and, at the low scale mz, 

M{mz) = MD{mo) + Ms{mz) , Ms{mz) = M^ad , (5.3) 

where Mrad is the matrix of radiative corrections. In general, Mj-ad can be written as 

Mrad = MsM(MSSM) + M^S , 

where M]^s is the contribution of non-standard corrections. In particular, we will analyze 
the corrections given by eq. (|4.8|), assuming ky^ = 0. In this case, Mrad can be written as 

Mrad = XiXMD + MdX^) , (5.4) 



where 



(5.5) 



The non-standard (flavor changing) couphngs Y affect also the evolution of the charged 
lepton Yukawa coupling matrix Yi, inducing non-zero off-diagonal entries. As a conse- 
quence, flavor basis should be redefined at the scale mz through a rotation Ui of charged 
leptons. Due to the strong hierarchy of charged lepton masses, the mixing angles gen- 
erated in the charged lepton sector are small, being of the order of radiative corrections 
themselves: Ui = 1 + C/rad- Therefore, the neutrino mass matrix in flavor basis is given by: 

M^' = Ui {Md + M,^)U{ ^Md + Mrad + U.^^Md + MdU^^ , (5.6) 

where we have neglected terms quadratic in radiative corrections. While the neutrino 
masses at mz can be extracted from M = Mjj + Mrad, the mixing angles are influenced 
by Ui rotation and all the terms in eq. (|5.6| ) should be considered. 

The term (C/rad-^D + -^-D^rad) describing the effect of charged leptons, 

has the same structure as Mrad in eq. ( |5.4| ). Therefore, the effect of rotation to flavor basis 
amounts to a redefinition: 

X ^ X + ^ , (5.7) 
A 

wheqre X and A are defined in eq. (^^- In the following, we will extract some results 
independent from the form of X. 

Let us comment on our assumption of exact zeros at mg (see eq. ( ^.2[ )). Zero val- 
ues of matrix elements can be a consequence of certain symmetry at ruQ. The radiative 
contributions to these elements are generated by interactions which break this symmetry. 
In general, the symmetry breaking leads to finite contributions already at mo- On the 
other hand, if some particles producing radiative corrections are light {nix *C mo), the 
largest contribution to zero elements is given by the leading logarithms and can be com- 
puted by RGE methods in the context of the effective theory. We are discussing here these 
logarithmic contributions. 

In what follows we will consider corrections to different dominant structures Mo- 

5.1 Unit matrix 

Let us consider the matrix with dominant structure proportional to the unit matrix: 

Mil = mdiag(l,l,l) =ml. (5.8) 

At mz, the approximate degeneracy of neutrino masses (mi ~ m2 ~ ma) is broken by 
Am^tjjj and Am^^j. Let us define 

„_1 rn.^ Aml^ ^ Amj, 

m2 2 m| m2 2m| 

where the sign of r] is + for inverted ordering (ms < m2) and — for normal ordering 
(ms > m2). Taking m^ 0.3eV (a value allowed by all existing upper bounds [§, |3|, ^]), 
one finds e <C 7? < 5 • 10"^ This means that deviations from M{mz) = ml can be smaller 
than 1%. 
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For simplicity, at mz we assume zero Majorana phases, sin^is = and maximal 
atmospheric mixing (^23 = Then, the phenomenological mass matrix can be written 

as \M: 



mz) 



1 - ecfa eci2Si2/\/2 -eci2Si2/^ 
m2 I ... 1 - (r/ + es?2)/2 (-r? + esf2)/2 



(5.9) 



... ... !-(?? + 6^^2)72^ 

where C12 = cos0i2 and S12 = sin 6*12. 

Let us consider, first, radiative corrections in the case of SM (MSSM). Substituting 
the matrix ( |5.8D in eq. ( |3.6| ), it is evident that off-diagonal elements remain zero and, 
therefore, no mixing is produced. Standard corrections lead to a split among diagonal 
elements. Using eq. ( p.6[ ), one finds that mass squared differences in the range required by 
phenomenology can be obtained at mz in the case of MSSM with large tan /?, however one 
gets the following prediction: 



m 
j 

m. 



3.5 • 10" 



This ratio is about ten times smaller than the best fit value ~ 2.5 • 10~^. 

Let us consider now radiative corrections which originate from new scalars and fermions 
(section ^T| ). Substituting the matrix ( ^. in eq. (^.4), we get corrections to matrix 
elements of the following form: 



-^rad ~ 2 m A 



Lee ReXe^ ReXe 



V 



(5.10) 



where X and A are defined in eq. ( |5.5| ). Since the matrix X is hermitian, the corrections 
to all matrix elements are real, so that no CP violation is induced in this case. 

The corrections to the different elements are related if a specific form of the matrix X 
is given. We consider an interaction like those in eqs. (|4.4| )-(4.6). Then, X^js 
defining 

Aq, = 2X\Ya\'^ , Cai3 = cos(arg - arg Yg) , a 



YaY^ and. 



we get 



M{mz) = Md + Mrad 



e,n,T, 



(5.11) 



m 



1 + A 



V 



(5.12) 



fj, \J X/l^TCfXT 

1 + Xr J 

The small parameters Ac^ are real and have all the same sign, which is opposite to the sign 
of ky^ (see eq. ( p.5D ). The parameters Cap can vary between —1 and 1, but only two of 
them are independent. 

Before comparing the matrix ( 5.12| ) with the phenomenological mass matrix ( |5.9D , one 
needs to perform an additional rotation Ui. If Ui is real, substituting Md oc 1 in eq. ( ^^ ) 
we get 

M^^ =Md + M,ad + 0(C/radM,ad) • 



In this case, we can neglect the rotation Ui. 
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Let us show that the matrix ( |5.12| ) can reproduce the matrix (|5.9| ). For simphcity, we 
take m = m2- Then, from the condition Md + M^ad = M'p^{Mz), we find: 



-eci2 



-2(^ + ^*12) 



1 
'2 



^1, 



12 



S12 ■ 



(5.13) 



Since Ag is negative, also A^ and A,- are negative. This corresponds to inverted mass 
spectrum (77 > 0). To satisfy eq. ( ^.13| ), one needs a positive ky^ . 

If m is not equal to m2, other solutions are possible. It turns out that the matrices (| 
and ( 5.12 ) can be equal only if m is one of the three eigenvalues. Requiring m = mi 



(1 — e)m2, one gets 

Ap ~ €'Si2 ) A,, 



A. 



1 
2' 



-efi 




V 



C12 • 



(5.14) 



This corresponds to normal mass spectrum. Requiring m = 



Ae ~ 



(1 — r/)m2, one gets 

eCi2Si2 



--en 



(5.15) 



This can correspond to both normal and inverted mass spectrum. 
5.2 Dominant M^e and M^r 

Let us consider the hierarchical matrix with dominant block given by 



Md 



m 



/l 

0-1 
(^0 -1 



For this matrix ^23 = 7r/4, ^13 = and the eigenvalues equal (m,m,- 



-m 



(5.16) 



Substituting 



the matrix (5.16|) in eq. (|3.6|), one sees that SM (MSSM) radiative corrections generate 



a solar mass squared difference: Am^^j ~ 2m?'er\k\. However, atmospheric mass squared 
difference and solar mixing angle remain zero. 

We assume that the relative phases among eigenvalues at mz are the same as at mo, 
that is (1, 1, —1) and also the values of ^23 and ^13 remain 7r/4 and 0, respectively. Then, 
the phenomenological mass matrix is real and can be written as pO|: 



( I- ecl^ eci2Si2/\/2 



mz) 



m2 



(ry 



es 



12 



)/2 



-eci2Si2/\/2 
■l + (^ + es?2)/2 



(5.17) 



v 



(ry 



12 



)/2 



where e, ry, C12 and S12 are defined as in section 5.1 
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Let us study non-standard radiative corrections. Substituting the matrix ( 5.16| ) in 
eq. ( ^.41 ), we get 



(5.18) 



1 + Ai As -A^ \ 

M(mz) = Md + M,ad = m I ... A3 -I + A4 

^ J 

where, using the definitions in eq. (|5.5|), 



fie 



-2 AX 



A4 



T 



To reproduce the phenomenological matrix ( |5.17] ) with the matrix ( |5.18| ), one needs to take 
A2 and A3 real (Ai and A4 are real by definition, since X is hermitian). Then, assuming for 
simplicity m = m2, the equation M{mz) = M'^^ijnz) leads to the following relations: 



tan 2^12 



2y2A2 



A3 — Ai 



A4 



r? = A3 + A4 



A4 — Ai — A3 



(5.19) 



YoYp, where are non-standard Yukawa 



Let us consider the specific form X^p 
couplings. In general, before a comparison with eq. ( 5.17 ), this form of X should be 
modified as in eq. (^.7|). Assuming that the effect of this redefinition is very small or can 
be reabsorbed in a redefinition of Ya, one can express the equality M{mz) = MP^{mz) 
in terms of the parameters defined in eq. ( ^.ll]) . In particular, the equality of matrices is 
realized for m = and the same values of parameters as in eq. (5.15). 



Non-standard radiative corrections to the matrix ( 5.16| ) has been discussed also in [p6|| , 
in connection with a non-abelian discrete symmetry that leads to the matrix ( [5. 161) at high 
energy. 

5.3 Dominant Mee and ^ur-block 

Let us assume that the matrix with dominant block at the scale niQ is given by 



m 



/l 

1/2 -1/2 
^0 -1/2 1/2 



(5.20) 



This matrix leads to 823 = 7r/4, ^13 = and the eigenvalues equal (m, m, 0). It corresponds 
to inverted mass spectrum. The SM (MSSM) radiative corrections do not generate 1 — 2 
mixing because the zero elements Mg^ and Mgr are not modified by these corrections. 

Let us consider the phenomenological mass matrix with the dominant block struc- 
ture ( ^.20 ). We assume, for simplicity, 1713 = arg(mi/m2) = ^13 = and 623 = it/A. Then, 



the matrix can be written as |20|: 



/ 1 - ecf 2 ecuSu/V^ -ecusu/V^ 



mz 



1712 



(1 



es 



\ 



i2)/2 -(1 

(1- 



12 



es 



12 



)/2 
)/2 



(5.21) 



where m2 = yA 



0.05eV and e Am^^^/ {2Am[ 



atm) 



0.01. 



-14- 



Let us study non-standard radiative corrections. Substituting the matrix ( 5.2C| ) in 
eq. ( ^.41 ) and requiring Mg^ = —M^t and M^^ = Mrr (necessary to reproduce the ma- 
trix (|5.21|) ), we get 



M{mz) = Md + Mrad = m 



where, using the definitions in eq. (|5.5|) , 

Ai = 2\Xee , A2 = 2AReXe^ , 



/l + Ai A2 -A2 

. 1/2 + A3 -1/2 -A3 
1/2 + A3 



(5.22) 



As = A(X 



Then, assuming for simphcity m = m2, the equation M{mz) = MP^{mz) leads to the 
following relations: 

2A3 



A2 — V -^1^3 



tan di2 



Ai 



-Ai - 2A3 . 



(5.23) 



If the matrix X has the specific form Xaf^ = YaY^ , one can express the equahty 
M{mz) = MP^{mz) in terms of the parameters defined in eq. (5.11). For m = m2, we get: 



A. 



-ec?2 



A/t — Al- 



es 



12 



-1, 



5.4 Dominant /ir-block 

Let us take, at the scale mo, the matrix with dominant /xr-block [p4 



Md = m 



1. 




(5.24) 



It corresponds to the case of normal hierarchical mass spectrum (7713 ^ m2 ^ mi). 

The mass matrix required by phenomenology can be written as an expansion over the 



small parameters cos 26*23, sin6'i3 and r = ^/Amg^j/Am^^j^ [19|: 
MP'^{mz 



2 



/ 
... 1 1 I + cos 2^23 



1 



+ sin 6*136' 



iS 



f V2 \/2\ 
... 

V / 



/ 
... -10 I + 

V 1 

/ 2sf2 \/2si2Ci2 -\/2si2Cl2 



(5.25) 



re 



2icr 



'12 



V 



"12 



-12 



where (5 is the CP violating Dirac phase and a is the CP violating relative Major ana phase 
between m2 and m^. 

Let us discuss, first, the effect of SM (MSSM) radiative corrections to the matrix ( 5.24| ). 
Using eq. ( |3.6|) , one sees that e-row elements remain zero along the RGE running. It is 
a consequence of eq. (^^). In contrast, SM (MSSM) corrections to ;ur-block elements are 
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present. These corrections induce a small deviation from ^23 = 7r/4, which can be easily 
computed using eq. (^): 



cos 29 



23 



{e^-er)k + 0{e^)^ 



The first neutrino remains massless and unmixed. Also ra2 remains zero because the deter- 
minant of the /xr-block is zero even after the inclusion of radiative corrections. Therefore, 
the solar mass difference and mixing angle are not generated radiatively. 

Let us consider, now, the effect of non-standard radiative corrections. Using the dom- 
inant matrix (5.24), the matrix of corrections, given by eq. (^.4|), is: 



rad 



m 



/ Ai Ai \ 
. . . 2A2 A2 + A3 

V 2A3 I 



(5.26) 



where, using the definitions in eq. (|5.5|) , 

Al = A(Xe^ + Xer) , A2 = A(X^^ 



\{Xr^, + X 

T 



The eq. (5^) implies that a given element Map receives corrections proportional only to 
elements in the a and fi rows of M/j. Since the first row and column in eq. ( 5.24 ) are 
zero, the element Mge remains zero. Corrections to M^,^ and Mer are equal. The matrix 
M = Md + Mrad has a zero eigenvalue, so that strong normal hierarchy is preserved by 
radiative corrections. 

Comparing the matrix of radiative corrections ( ^.26 ) with the phenomenological ma- 
trix ( 5.25| ), we find that the atmospheric angle gets a deviation from 7r/4 and the angle ^13 
becomes non-zero: 



cos 26*23 ~ A3 - A2 , sin 6*13 Ai 



V2 • 



(5.27) 



A problem appears with the generation of solar parameters, because the structure of the 
term proportional to r in eq. ( ^.25] ) (let us call it M^°^) differs from the structure ( 5.26| ) of 
radiative corrections: 



M: 



sol 



-M: 



sol 



sol 

/iT 



-M; 



sol 



-M: 



sol 



whereas 



Mlf = , 



M. 



rad 



M: 



rad 



M. 



rad 



rad\ 

TT ) 



In fact, computing the eigenvalues of M = Mo + Mrad, one finds that r is of order A?: 



7712 

772-3 



^"^atm 



A? + (A2 - A3)V2 



(5.28) 



To satisfy the lower bound (99% C.L.) on r, which is ~ 0.1 (LMA), one needs A, > 0.2H-0.3. 
In other words, LMA solar mass difference can be generated only if radiative corrections 
are ~ 10 times larger than in the MSSM with large tan (3. 
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The equahties Mge = and Mg^ = Mer(= mXi) in eq. ( ^.26| ) lead to the following 
general relations between observables [19|: 

2 = ^ ^ ^ (1 - sin 2023) ) tan 6*13 = si2\/r, sin(5 = 0, (5-29) 

The first equality in ( 5.29| ) shows that the LMA solar mixing angle can be obtained only 
if atmospheric mixing deviates significantly from maximal value and if r is substantially 
smaller than the present best fit value. The equality can be satisfied taking 99% C.L. 
allowed intervals for 9i2, 023 and r. The other equalities in ( ^.291 ) show that the radia- 
tive corrections ( |5.26| ) give S13 close to the present upper bound ~ 0.2 and lead to CP 
conservation in oscillations and to opposite CP parity between m2 and 771,3. 

We have shown that, using the form ( 5.26 ) of corrections, the predictions for ^rn'^^^ 
and 9i2 are too small with respect to phenomenological values. However, the corrections 
related to the second term in square brackets of eqs. ([4.^), ( 4.11 ) have qualitatively different 
features and can give better predictions for the solar parameters. In particular, a non-zero 
Mee can be generated by these corrections. 

For example, let us consider the interaction in eq. (|4.2|). The /3-function coefficients 



are, in this case, ky^ = 1/2 and fey^ = 1. The matrix of couplings Y is antisymmetric, 
therefore the matrix of corrections AM^ (see eq. ( [4.8D ) depends on three independent 
couplings only: 1^^, Ygr and Y^r- We want to compare M{mz) = Md + AM^ with 
the phenomenological mass matrix ( 5.25| ) (we neglect, for simplicity, the charged lepton 
rotation Ui). The equality of the two matrices can be realized for sin0i3 = cos 2^23 = 
and 

^-^ + ^' Ay,^ - Ay,, - -y , Ay^,-^^^— -y, (5.30) 

where A is defined in eq. (^). 

Notice that, while the form (|5.26| ) of corrections predicts r of order A?, here r is of 
order Aj. Therefore, ~ 10% radiative corrections are enough to generate radiatively Am^^j. 
Such corrections can be produced if three generations of new particles (e.g., three scalar 
singlets) are introduced. As follows from eq. ( 5.30| ), also the solar mixing angle 9i2 can be 
generated radiatively (in the LMA allowed range). 



6. Discussion and conclusions 

The structure of the Majorana neutrino mass matrix can be reconstructed, using experi- 
mental data, at the electroweak scale mz- This structure changes with RGE running to 
the high energy scale mo- We have analyzed the features of the running between the two 
scales. 

The SM (MSSM) radiative corrections do not modify the matrix structure. In flavor 
basis, the value of each matrix element at mz is proportional to the one at rriQ. Moreover, 
the corrections to this value cannot be larger than few percents. Therefore, both the 
dominant structure of the mass matrix and the small matrix elements are not modified 
significantly between m-o and mz- Zero elements remain zero. 
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At the same time SM and MSSM corrections can change significantly observables: 
corrections can enhance or suppress mixing, modify strongly Am^ or even generate mass 
split. Substantial change of observables occurs for the quasi-degenerate spectrum, with 
common scale of neutrino mass m > 0.1 eV. 

We have studied the radiative effects induced by new particles and interactions at a 
scale nix, with mz < mx < niQ. These non-standard (fiavor changing) corrections lead to 
coupled RGE's of different matrix elements. As a consequence, small tree-level elements get 
corrections proportional to the large matrix elements. We have considered non-standard 
corrections induced by new scalar bosons, new fermions and new gauge bosons. 

In all cases, the dominant structure of the mass matrix remains the same between 
mz and mo- Therefore, if the matrix structure at mo is determined by some symmetry, 
this symmetric structure can be identified from the experimental data at mz- However, 
the values of small elements can be strongly modified. We show that small elements of 
hierarchical matrices can be zero at the scale mo and receive non-zero contributions from 
radiative corrections. At the high mass scale, only the dominant block elements can be 
non-zero. 

In the case of exactly degenerate neutrino masses, small (~ (0.1 -r- 1)%) corrections to 
zero elements can generate large mixing angles and mass squared differences in the range 
required by phenomenology, both for solar and atmospheric neutrinos. We have shown, in 
particular, that the unit matrix can be the exact form of the neutrino mass matrix at mo. 

In the case of inverted hierarchy, the structure with zero e/i and er elements at mo can 
lead to correct predictions for low energy solar parameters, if 1% non-standard corrections 
are present. 

In the case of normal hierarchical neutrino masses, we have studied the matrix with 
dominant /ir-block. Solar mass squared difference and mixing angle can get large renormal- 
ization effects. To generate a mass difference in the LMA region, one needs 10% radiative 
corrections. 
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